Excitation of a Dipole Topological Mode in a Strongly Coupled Two-Component 

Bose-Einstein Condensate 
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Two internal hyperfine states of a Bose-Einstein condensate in a dilute magnetically trapped gas 
of 87 Rb atoms are strongly coupled by an external field that drives Rabi oscillations between the 
internal states. Due to their different magnetic moments and the force of gravity, the trapping 
potentials for the two states are offset along the vertical axis, so that the dynamics of the internal 
and external degrees of freedom are inseparable. The rapid cycling between internal atomic states 
in the displaced traps results in an adiabatic transfer of population from the condensate ground 
state to its first antisymmetric topological mode. This has a pronounced effect on the internal Rabi 
oscillations, modulating the fringe visibility in a manner reminiscent of collapses and revivals. We 
present a detailed theoretical description based on zero-temperature mean-field theory. 

PACS numbers(s): 03.75.Fi, 05.30.Jp, 32.80.Pj, 74.50.+r 



I. INTRODUCTION 

An intriguing aspect of Bose-Einstein condensation 
(BEC) in a dilute atomic gas is that the internal atomic 
state of the condensate can be manipulated to produce 
quite novel systems. A number of interesting experi- 
ments have produced uncoupled multi-component con- 
densates, in which two or more internal states of the con- 
densate exist together in a magnetic or optical trap jij-fl . 
These experimental studies, along with their theoretical 
counterparts, have investigated various topics such as the 
ground state of the system , the elementary excita- 

tions ||12| — |lST , and the nonlinear dynamics of component 



separation |H 

Many fascinating properties can be studied by apply- 
ing an external electromagnetic field that coherently cou- 
ples the internal atomic states of the condensate S . 
In the experiment described in [ pof , the relative phase be- 
tween two hyperfine components was measured using a 
technique based on Ramsey's method of separated oscil- 



lating fields 1 29 1 and these results have motivated further 
theoretical investigation [^9|,|lJ. In this system, several 
key parameters can be varied over a wide range of val- 
ues, such as the coupling-field intensity and frequency, 
the confining potentials, the total number of atoms, and 
the temperature, making this a very rich system to ex- 
plore. 

Several theoretical papers ]22]-|2~4[] have investigated the 
weak coupling limit of this system, where the intensity 
of the external driving field is very low but is turned 
on for a time long compared to the period of oscilla- 
tion in the magnetic trap. A clear analogy exists be- 
tween this system and the Josephson junction f30|| . In the 
Josephson junction, identical particles in spatially sepa- 
rated condensates are coupled via the tunneling mecha- 
nism j39|l . In this two-component system, however, 
two distinct internal states of the condensate are coupled 
by an applied field. By adjusting the magnetic fields that 



confine the atoms, the degree of spatial overlap between 
the two components can be controlled. 

Recently, there have been both experimental and 
theoretical ||(J studies on the dressed states of a driven 
two-component condensate, drawing an analogy to the 
dressed states of a driven two-level atom in quantum op- 
tics. Due to the interplay between the internal and ex- 
ternal degrees of freedom, the condensate dressed states 
have spatial structure that depends on the trap param- 
eters, the mean-field interaction, and the frequency and 
intensity of the driving field p6[ . In the experiment re- 
ported in 1 25 1, the dressed states were created via an 
adiabatic passage by sweeping the detuning. 

In this paper, we focus on the limit of a very strong and 
sustained coupling between hyperfine states, which is the 
situation achieved experimentally in p5[ . In that exper- 
iment, a BEC of about 8 x 10 5 atoms was produced in 
the \F = 1, Mp = —1) hyperfine state of 87 Rb, at a tem- 
perature close to zero, T » 0. The atoms were confined 
in a time-averaged, orbiting potential (TOP) magnetic 
trap by a harmonic potential with axial symmetry along 
the vertical axis. An external field was then applied that 
coupled the |1, —1) state to the |2, 1) hyperfine state via 
a two-photon transition. The Rabi frequency was five 
to ten times larger than the vertical trap frequency and 
the detuning could be adjusted arbitrarily. Due to their 
different magnetic moments and the force of gravity, the 
two hyperfine states sit in shifted traps offset along the 
vertical axis jl0| . The degree of separation could be con- 
trolled by adjusting the magnetic trapping fields. 

The subsequent behavior of the system described in 
p5[ was quite unexpected: after the coupling field was 
turned on, the Rabi oscillations between the states ap- 
peared to collapse and revive on a time scale which was 
long compared to the Rabi period of 3 ms. An example of 
this behavior taken from p5| is shown in Fig. 1. It was 
observed that the period of this modulation increased 
with decreasing detuning. The behavior of the system 
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also depended critically on the separation between the 
traps for each state. As the separation was taken to 
zero, the effect went away. The evolution of the density 
of each component was also very interesting. Each com- 
ponent cycled between a density profile with one peak 
and a profile with two peaks. The two-peaked structure 
was most clearly visible around the collapse time, which 
is 1 1=3 20 ms for the case shown in Figure 1. 
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FIG. 1. This plot shows the modulation of the fractional 
population in the (1,-1) state. The top line is experimen- 
tal data |^| while the bottom line is the result of a nu- 
merical calculation of the three-dimensional, two-component 
Gross-Pitaevskii equation Eq. ([!]) (below). The coupling 
strength and detuning were chosen for the calculation to be 

= 350 Hz and 8 = —188 Hz, respectively. 

A numerical calculation of the three-dimensional, cou- 
pled Gross-Pitaevskii (GP) equations Eq. (|l]) (below) de- 
scribing the system in the zero temperature limit agrees, 
at least qualitatively, with the outcome of the experi- 
ment, as shown in Figure 1. Of course, the agreement 
between the numerical integration of Eq. ([!]) and the lab- 
oratory data does not in itself provide an intuitive expla- 
nation of the underlying physical mechanism responsi- 
ble for the collapse-revival behavior. In this paper, we 
present a detailed analysis of this problem, and arrive at 
a rather simple model that explains the major features 
of the system's behavior. 

Before presenting the details of our analysis, it is use- 
ful to first give an overview of the results. There are two 
main concepts that play key roles in obtaining an intu- 
itive understanding of this problem. First, there is a clear 
separation of time scales: the period of Rabi oscillations 
between internal states is much shorter than the period 
of the trap. That is, the internal dynamics occur on a 
much shorter time scale than the motional dynamics of 
the system. It is therefore useful to go to a frame rotating 
at the effective Rabi frequency. In this rotating frame, we 
show that there exists a weak coupling between the low 
lying motional states which is proportional to the offset 



between the two traps. This weak coupling has the effect 
of modulating the amplitude of the fast Rabi oscillations 
in the lab frame. 

The second key point is to understand exactly which 
motional states are excited. They are not the linear re- 
sponse collective excitations (normal modes) that have 
been studied frequently in the BEC literature [pd|-|43|. 
Instead, they are many-particle topological modes de- 
termined by the self-consistent solutions to the two- 
component GP equations Eq. (|l4|). The well known vor- 
tex mode p4|-[i7[ is one example of such an excitation 
in which phase continuity requires quantized circulation 
around a vortex core. The related excitation which plays 
a key role in this paper does not exhibit circulation but 
has a node in the wavefunction amplitude and exhibits 
odd-parity behavior characteristic of such a dipole mode. 
For a single component in the limit of the uniform gas, 
the exact solution of this mode is known as a dark soli- 



ton 1 39 1 . In that case the scale of the density perturbation 
around the node is the healing length and is determined 
by a balance of kinetic and mean-field interaction en- 
ergies. In the problem we consider here, however, it is 
necessary to account for the mean field of the the remain- 
ing population in the condensate ground state. Conse- 
quently the two modes — the ground state and the dipole 
mode — are inextricably linked and must be determined 
self consistently. 

We first present a detailed theoretical analysis in Sec- 
tion H. After making several reasonable approximations, 
we arrive in Section II E at the two- mode model - a sim- 



plified description that encapsulates the essential prop- 
erties of the system. In Section [II we present results of 



numerical calculations that illustrate the behavior of the 
system and we compare our model with the exact nu- 
merical solution of the coupled GP equations. We finally 
summarize our work in Section ^ and suggest further 
studies based on our understanding of this phenomenon. 



II. THEORETICAL DESCRIPTION 

The following theoretical development was motivated 
by the experiment described in p5|| . Therefore, we have 
not tried to keep our calculations general, but instead 
have made several assumptions based on that particu- 
lar situation. However, our approach could easily be ex- 
tended to treat a broader class of systems. We give a brief 
discussion in the conclusion of the paper about possible 
extensions of this work to other interesting systems. 

We begin this section by writing down the coupled 
mean-field equations, valid for zero temperature, that de- 



scribe this driven, two-component BEC. In Section [IB 



we rewrite the mean-field equation in a direct-product 
representation that clearly separates out the external and 
internal dynamics. We then go to a frame rotating at the 
effective Rabi frequency in Section II C in order to focus 
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on the slower motional dynamics of the system. After 



making some approximations in Section [ID, we finally 
arrive at the main result of our study in Section [I E , the 
two- mode model. 



A. Coupled Mean-field Equations 

A mean-field description of this many-body system 
that includes the atom-field interaction has been devel- 
oped, which generalizes the standard Gross-Pitaevskii 
(GP) equation to treat systems with internal state cou- 
pling H|,|2^,Q . The resulting time-dependent GP equa- 
tion describing the driven, two-component condensate is 
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H° + Hf lF + 5/2 
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fi/2 



(i) 



The Hamiltonians describe the evolution in the trap H® 
and the mean field interaction iT 2 MF for each component 



I 2 a 



(2) 



where 71 = — 1 and 72 = 1, and zq is the shift of each 
trap from the origin along the vertical axis. The factor 
a = lu z /lu p is the ratio of axial and radial trap frequen- 
cies. In the experiment reported in pf|, a > 1. The 
mean-field strength is characterized by Xij = 47ray/z s ho, 
which depends on the scattering length ay of the colli- 
sion. In general there will be three different values, one 
for each type of collision in this two-component gas: an, 
0-22, a-12- The detuning between the driving field and 
the hyperfine transition is given by S while SI denotes 
the strength of the coupling. We work in dimensionlcss 
units: time is in units of 1/uj z , energy is in terms of 
the trap level spacing %u) z , and position is in units of 
the harmonic oscillator length z sn o = V^/VriRb u z- The 
complex functions ipi (r, t) are the mean- field amplitudes 
of each component, where i = {1, 2}. They obey the nor- 
malization condition J (\ipi\ 2 + \ipj\ 2 )d 3 x = 1. The total 
population is N. 

The coupled mean-field equations Eq. (|l|) can be solved 
numerically using a finite-difference Crank- Nicholson al- 
gorithm ^9|. We show results of such calculations in 
Fig. 1 for the three-dimensional solution and in Sec- 
tor a one-dimensional version. However, in order 
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to gain a more intuitive understanding of the behavior 
shown in Fig. 1, we formulate a simplified description of 
the system in the following section. 



B. External ® Internal Representation 

The coupled mean-field equations Eq. (|l|) can be 
rewritten in a more illuminating form by making a clear 
separation of the external and internal degrees of free- 
dom. The system exists in a direct-product Etilbert space 
T~l = Hex ® Win, where TC ex is the infinite-dimensional 
Etilbert space describing the motional state of the sys- 
tem in the trap and 7ii n is the two-dimensional Etilbert 
space describing the spin of the system. A general opera- 
tor in Ti. can be written as a sum over the direct-product 
of operators from H cx and H m . We rewrite Eq. (|l]) in 
this representation as 

d f2 S 

* 5jj£llK*)} = [#0 ® i + i ® (-j °* + +h z <8> v z ]\i>(t)) 

(3) 

where {1, a x , & y , a z } are the standard Pauli spin matri- 
ces. The state of the system \ip(t)) in general has a 
nonzero projection on the internal states |1) and |2), 
represented by ipi(r,t) — (r\(i\ip(t)), where i = {1,2}. 
The position representations of Hq and H z are local, i.e. 
<r|£ |r') - H (r)S(r - r') and (r|^ z |r') = H z (r)5{v - 
r'), where Ho(r) and H z (r) are given by 

Ho(r) = -\V 2 + \[(ty + z 2 ] + {^{t)\P r ® X + m)) , 
2 2 a 

H z (r)=-z z+Uj(t)\P r ®\^(t)) . (4) 

The operator P r is the projector onto the position eigen- 
states P r = |r)(r|, and the matrix representations of A+ 
and A_ are given as 



j = NfX 1 + \ 12 
+ 2 1 A 2 + A 12 



A_ 



N 
~2 



Ai - A 12 
A12 — A 2 



(5) 



Note that the harmonic potential in Hq is centered at the 
origin. The mean-field interaction has been rewritten in 
terms of a part that acts identically on both components 



|Pr«)A4 



) 1, and a part that acts with the opposite 



sign on each state (ip\P r ® A_h/>) ® <j z . 

The first two terms in Eq. (g) separately describe the 
external and internal dynamics of the system, respec- 
tively. The third term in Eq. (^), however, couples the 
internal state evolution to the condensate dynamics in 
the trap and can lead to interesting behavior. If the 
term H z were identically zero, then the problem would 
be completely separable in terms of the external and in- 
ternal degrees-of-freedom. The term H z would be zero 
if the trap separation Zq were zero and if the scattering 
lengths were all exactly equal. In fact, for 87 Rb the three 
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scattering lengths are nearly degenerate, so the main ef- 
fect of H z comes from the term — zqZ, which is the differ- 
ence in the shifted traps. It causes there to be a spatially 
varying detuning across the condensate. 



C. Rotating Frame 

As previously stated, we are concentrating on the sit- 
uation where the coupling strength is large, so that the 
frequency of the Rabi oscillations f2 is significantly larger 
than the trap frequency v z . In this case, the internal spin 
dynamics and the motion of the condensate in the trap 
occur on two different time scales. Therefore, it is useful 
to go to a rotating frame that eliminates the second term 
in Eq. @ describing the fast Rabi oscillations between 
the two internal states. In the rotating frame, we will 
be able to understand more clearly how the third term 
in Eq. (|J), which couples the motional and spin dynam- 
ics of the condensate, effects the system on a time scale 
much longer than the period of Rabi oscillation. 

We go to the rotating frame, or interaction picture, by 
making a unitary transformation using the operator 



Uj(t) 



- e -»i®(4 



(6) 



This can be rewritten in the equivalent form 



Ui(t) = 1 <g> (cos(fi e ff/2t)i - 7— sm(n cS / 2 t)[Qa x + 5a z 

"eff 

(7) 



where f2 c ff = \/ SI 2 + <5 2 . The state vector l^Ht)) m the 
rotating frame is related to the state vector in the lab 
frame \ip{t)) by 

\i>w{t)) = u\m)) . (8) 

In the rotating frame, the system evolves according to 



i^ I \t)) = H^(t)\^ I \t)) 



(9) 



where H^ (t) is the interaction Hamiltonian 

H {I) (t) = H <8> i + H z <8> (a x (i)<7 x + a y (f)oy + a z (t)a z ) . 

(10) 

Note that Ho and H z are unaffected by the unitary trans- 
formation to the rotating frame. The time-varying coef- 
ficients a x (t), a y (t), and a z (t) are 



*x(<) = 



n 



[1 - cos(fi off i)] 



a y (t) = — — sin(ficffi) 

"off 



S 2 

az (*) = Q2~ + Q2~ COs(fiefft) 
eff " " 



^efit 



(11) 



D. Approximations 

We now make three simplifications in order to ex- 
tract out the dominant behavior of the system. We first 
note that the coefficients oti{t) given in Eq. (^Tj) oscil- 
late rapidly at the Rabi frequency. We expect the sys- 
tem in the rotating frame \ip^(t)) to evolve on a much 
slower time scale than the period of Rabi oscillation. We 
can utilize this fact in order to simplify the interaction 
Hamiltonian H^ (t) given in Eq. ([n]) by taking the av- 
erage values of the coefficients this is equivalent 
to coarse graining Eq (Q). The coefficients in Eq. ( p"l| ) 
become time independent and reduce to: a x = <5f2/f2^, 
ct y = 0, and a z = 5 2 /Og ff . 

We make the additional assumption that the system 
is being driven close to resonance, so that <5/fi c g << 1. 
We therefore set a z = since a z depends quadratically 
on this small parameter. 

Finally, we take advantage of the fact that the scat- 
tering lengths for 87 Rb are nearly degenerate, with the 
ratios between inter- and intra-species scattering lengths 
given by {a 2 : a 12 : 01} = {0.97 : 1 : 1.03} [|. 
This allows us to simplify the two mean-field terms ap- 
pearing in Eq. (Q). We first make the approximation 
A + ~ XN • 1 by assuming equal scattering lengths, or 
Ai w A2 ~ A12 — > A. We can also simplify the other 
term A_ by assuming that its predominant effect is to 
) , shift the levels slightly. Instead of neglecting it alto- 
gether, we simply replace it by a mean-field shift of the 
levels (i/j\P r <S> A_|i/)) — * Smf, where the shift is given 
by Smf = AA J n 2 (r,0)d 3 x/N. Here n(r, 0) is the total 
density at t = and AA = (Ai — A12) = (A12 — A 2 ). 
We can absorb it into the detuning 5 by defining an 
effective detuning 6' that includes this mean-field shift 
5 — > 5' = S + 5mf- 

After making the above approximations, we can now 
write the interaction Hamiltonian H W from Eq. ( |To| ) in 
a much simpler form 



H' 



i + H' 



(12) 



where the position representations of H' and H' are lo- 



cal, i.e. 



<r|W 



H^{r)S(r - r') and (r\H' z \r') 



H' z (r) 5(r — r'), where Hq(t) and H' z (r) are given by 



H (r) 



2 

-0z 



V 2 + ^(-) 2 + 2 2 ] + An(r,t) , 
a 



(13) 



The total density is n(r,t) = N (ip(t)\P r ® i|-0(*))> and 
(3 = zoSfl/ft 2 s . For the typical values of the param- 
eters in the experiment, £1 — 400 Hz, 5' = 100 Hz, 
zo/zsho = 0.1, this factor is rather small (3 « 0.02 in har- 
monic oscillator units. Note that H' still varies slowly in 
time through the nonlinear mean-field term, which de- 
pends on the density. We refer to this result Eq. (Ill) 
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as the coarse-grained, small detuning (CGSD) model to 
distinguish it from the two-mode model presented below, 
which makes further assumptions. 

We have managed to greatly simplify the description 
of the system by going to the rotating frame. The first 
term in Eq. (|l^) contains the kinetic energy, a harmonic 
potential centered at the origin, and a mean-field interac- 
tion term depending on the slowly varying density n(r, t). 
The second term in Eq. ( |l2| ) represents a very weak cou- 
pling between the two internal states |1) and |2), and 
between motional states |</> n ) and \4>m) via the dipole op- 
erator z. The states \4> n ) and \4> m ) are the instantaneous 
self-consistent eigenmodes of H' . In the next subsec- 
tion we present a model that assumes only two motional 
states are coupled, the self-consistent ground state |0o) 
and the self-consistent first-excited state \4>i), which has 
odd-parity along the z-axis. 



E. Two-mode model 

It is useful to define a basis of motional states with 
which to describe the system in the rotating frame. A 
natural choice is the set of instantaneous eigenstates of 
Hq, which satisfy 



i-^ 2 + k-) 2 + A + An(r,t))^(r) = e, &(r) 



2 LV a' 



<p i (r)c)> j (r)d 3 r = 5 i , j , (14) 



where the index i refers to all of the relevant quan- 
tum numbers that uniquely specify each eigenstate, i = 
n p , n^,}, given the cylindrical symmetry of the sys- 
tem. In general, many modes can be occupied and the 
state vector is written 

\^ i Ht)) = Y^[a(tMi)\i)+d i (tmm} , as) 



where <j>i(r) 
is then 



(r|0i) . The density appearing in Eq. ( M ) 



n(r,t) - N(\J2^(t)Mr)\ 2 + | ^>(*)<Mr)| 2 ) . (16) 

i i 

It is clear that the set of coupled eigenvalue equations 
given in Eq. ( |l4| ) is nonlinear and requires a numerical 
procedure that will converge upon the solution in a self- 
consistent manner. The eigenstates (j>i(r) and eigenener- 
gies ti depend on time implicitly through the coefficients 
Ci (t) and di (t) , however we do not show this time depen- 
dence in order to simplify the notation. We assume that 
the eigenbasis evolves slowly in time so that the adiabatic 
condition is satisfied [f50| . 

Based on the experiment reported in p5|] the initial 
motional state of the system is ^^(r,0) = 0o(r — zqz); 



the system is in the ground state of Hq, but displaced 
from the origin along the vertical axis by zq. This dis- 
placement is small compared to the width w z of the con- 
densate zq/w z Rj 0.01. We therefore approximate the 
initial state of the system as \ip^\t)) — |</> )|1). 

The system in the rotating frame evolves according 
to the Hamiltonian described by Eq. ( |l2"| ) and Eq. (|l3|). 
The term H' z a x couples the internal states |1) and 
1 2) via <r x . It also drives transitions between motional 
states via the dipole operator z. The dipole matrix ele- 
ment (z)ij = {4>i\z\<pj) is the largest between neighboring 
states and falls off quickly as \i — j\ increases. For a small 
coupling parameter /3, we expect the coupling to the first 
excited state \<fix) to dominate the other transitions, mak- 
ing the evolution of the system predominantly a two state 
evolution. We therefore make the approximation that the 
system occupies only two modes 



|^ (/) (i)> =co(t) |0o)|l> +rfi(t)|0 1 >|2) 



(17) 



where |</>o) is the ground state i — {0,0,0} and |0i) is 
the first excited state with odd parity along the z-axis 
i = {1,0,0}. 

If we substitute this ansatz into Eq. (||), using the 
Hamiltonian described by Eq. ( |l2| ) and Eq. (|l3|), we get 
the equation of motion for the coefficients Cq (t) and di (t) 



-0{z)oi 



-P{z)oi 



(18) 



where we have neglected the time rate-of-change of the 
slowly varying adiabatic eigenbasis. This coupled pair 
of equations must be solved numerically by updating the 
energies &j and the dipole matrix element (z)oi from solv- 
ing Eq. (14) at each time step. However, in order to 



see how the behavior depends on the various physical 
parameters, one can obtain a simple estimate of the so- 
lution by fixing and (z)oi to their initial values. In 
this case the solution of Eq. (|l8|) is trivial and is given 
by co(t) = cos(rioi/2i) — i(Ae i/£l i) sin(O i/2 1) and 
di(t) = —i(2[3(z)/rioi)sm(Qoi/2t), where Aeoi = t\ — eo 
and Oqi = -\/4/3 2 (z) 2 + Ae^. In the rotating frame, the 
system oscillates between the two states at a frequency 
of Ooi: which is much slower than the effective Rabi fre- 
quency fieff- 

The oscillation frequency l^oi increases with increasing 
detuning 8' and increasing trap separation zq through the 
coupling parameter (3. The amplitude of oscillation de- 
pends on the energy spacing between modes Aeoi • Based 
on numerical calculations, we have found that this effect 
is enhanced by the mean-field interaction because Aeoi 
decreases with increasing population N. Also, the dipole 
matrix element (z) increases with increasing N, since the 
width of the condensate increases with increasing popu- 
lation. 

The solution in the lab frame can be obtained by ap- 
plying Ui{t) from Eq (g) to in Eq. @ to yield 
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m)) = ( ai (t) CQ (t) i^o) 

+ M*)c&(f) |<M 



■a2(t)di(t)|0i))|l) 
-ol(t)di(*)]0i»|2> 



(19) 



where cvi(t) = cos(O e ff/2i) — i(S'/Cl e s) sin(f2 e ff/2i) and 
"2(0 = — i(fi/fi e ff) sin(r2 e ff/2i). Eq. ( |T^ ) is the main 
result of our study, with which we can explain the essen- 
tial properties of the system. During the first few Rabi 
cycles t ss 1 / £l e ff , the coefficient di(t) 0, so that the so- 
lution for short times is \ip(t)) = (ai(i)|l)+a 2 (t)|2)) |0 O >- 
That is, for short times, the internal and external degrees 
of freedom appear to be decoupled and the system sim- 
ply oscillates rapidly between internal states. However, 
for longer times, the coefficient d\(t) grows in magni- 
tude as Co(t) correspondingly decreases. This results in 
a modulation of the Rabi oscillations. Furthermore, a 
two-peaked structure in the density appears, associated 
with the first-excited state \<j>i). 



III. RESULTS 

The main goal of this section is to illustrate the behav- 
ior of the system by showing results of numerical calcu- 
lations. For this purpose, it is useful to treat the system 
in only one dimension — along the vertical axis p^| . We 
also assume equal scattering lengths throughout this sec- 
tion, so that Smf = 0. Values of most of the physical 
parameters are given in Table 1. Values of the remain- 
ing parameters are stated for each case considered in the 
text. 



A. Understanding the dual dynamics 

In Figure 2 we plot the fractional population of state 
|1), given by Ni(t) = J | (z\ (l\ip(t)) \ 2 dz, for the case of 
fl = 700 Hz and S = 100 Hz. This is a numerical solution 
of Eq. (0) . The population is cycling rapidly at the effec- 
tive Rabi frequency Q c g = 707 Hz, while simultaneously 
being modulated at a much lower frequency of about 11 
Hz. 

In order to visualize how the spin and motional dy- 
namics become entangled over a time long compared to 
the Rabi period, we show snapshots of the density of each 
state in Figure 3. Three different sets of snapshots are 
shown, corresponding to the three circled numbers in Fig- 
ure 2. A full Rabi cycle is shown for each set. The first 
set begins at t = with all of the atoms in the 1 1) internal 
state and in the mean-field ground state of the trap \cj>o). 
During this first Rabi cycle, the shape of the density pro- 
file for each internal state does not change much — only 
the height changes. That is, the motional state remains 
the ground state while population cycles rapidly between 
internal states, as discussed below Eq. (|l9|). 




100 



40 60 
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FIG. 2. This plot shows an example of the modulation of 
the Rabi oscillations. The fractional population in state |1) 
is plotted as a function of time, obtained from a numerical 
solution of the one dimensional version of Eq. (Q). The values 
of the various parameters are given in the text. In Figure 3, 
the densities for both states are shown for three different Rabi 
cycles designated by the circled numbers in this plot. 








position (arbitrary units) 

FIG. 3. This plot shows snapshots of the density of each 
state for three different Rabi cycles corresponding to the three 
circled numbers in Figure 2. The solid line is the density of 
the |1) state, while the dashed line is that of the |2) state. 
Each snapshot within a set is numbered in sequential order. 
The first set starts at t = ms, and runs for a full Rabi cycle 
1.41 ms. The second and third sets begin at t — 45.2 ms 
and t = 90.3 ms, respectively. The time increment between 
snapshots is At = 0.28 ms for all three sets. 
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FIG. 4. The fractional population of the |1) state in the 
rotating frame is shown. The solid line is the solution given 
by the CGSD model, while the dot-dashed line corresponds to 
the solution of the two-mode model. If the two-mode model 
is extended to include coupling to the first even-parity excited 
mode, then we get better agreement to the CGSD model, as 
shown by the dashed line. 
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FIG. 5. The top strip of this plot shows snapshots of the 
density of each state corresponding to the solution of the 
CGSD model given by the solid line in Figure 4. The bottom 
strip shows the corresponding two self-consistent eigenmodcs 
given by the solution of Eq. (|14j). The times of each snapshot 
are shown in the region between the two strips, in units of 
milliseconds. The solid line corresponds to the density of the 
|1) state, while the dashed is that of the 2) state. 

The second set of snapshots in Figure 3 is taken at 
around t = 45 ms, which is halfway through the modula- 
tion. The density profiles for each spin state cycle rapidly 
between a single-peaked and a double-peaked structure. 
For example, in the first snapshot, the |1) state is in 
the single-peaked structure, while the |2) state is in the 
double-peaked structure, but halfway through the Rabi 
cycle the situation is reversed, as shown in the third and 
fourth snapshots. Finally, at about t = 90 ms when the 



amplitude of the Rabi oscillations has revived, the third 
set shows that the motional and spin degrees of freedom 
appear to be decoupled again, with the density profile of 
each spin state appearing as it did during the first Rabi 
cycle. 

As outlined in Section II, this peculiar behavior is most 
easily understood by going to the rotating frame. In 
Figure 4, we plot the fractional population in the 1 1) state 
in the rotating frame N [^{t) = j \(z\(l\^(t))\ 2 dz. The 
solid line corresponds to the CGSD model presented in 
Section IID. In the rotating frame, population is slowly 
transferred out of the |1) state due to the coupling from 
H' z ® a x in Eq. ([l|). 

In the rotating frame, the system is being excited out 
of the ground state |<^n) due to the dipole coupling H' z . 
This can be seen in the top strip of snapshots in Figure 
5, where the density of each spin state in the rotating 
frame is shown, corresponding to the solid line in Figure 
4. Initially, all of the atoms are in the |1) internal state 
and the mean-field ground state of the trap \4>o)- Due to 
the dipole coupling, population is transferred out of the 
ground state. 

The strongest coupling is between the ground \<fio) and 
the first excited \4>i) modes. These eigenmodes are shown 
in the bottom strip of Figure 5. They evolve slowly in 
time as the coefficients cn(i) and d\(t) change. For exam- 
ple, initially the ground state is just the Thomas-Fermi- 
like ground state, since all of the population is in that 
state. However, at t = 45 ms, about one-third of the 
population is in the first excited mode, which pinches 
the ground state due to the mean-field interaction term 
in Eq. (|ll|). That is why the self-consistent ground state 
at t = 45 ms is narrower than at t = 0. 

It is clear from Figure 4 that the low-frequency mod- 
ulation of the rapid Rabi oscillations in the lab frame 
is just the frequency of oscillation in the rotating frame 
between |0n)|l) and |(/>i)|2). This is reflected in the two- 
mode solution given by Eq. ( ft9| ) , which also helps explain 
the peculiar behavior of the densities shown in Figure 3. 
In the lab frame the system is cycling rapidly between the 
two modes shown in Figure 5. The initial values of the en- 
ergies are eg — 13.6 hu> z and e\ = 13.7 fko z , which makes 
Ani =0.1 hw z . This small energy splitting is due to the 
effect of the mean field, since in the limit TV — > 1 these 
energies move apart by a factor of ten, which greatly re- 
duces the coupling between the modes and thus greatly 
reduces the modulation effect. 

If we make the two- mode ansatz and solve Eq. (H), 
we get the dot-dashed line in Figure 4. The discrepancy 
from the solid line arises due to a weak coupling between 
the first |</>i) and second excited modes. If we extend 
our two-state model to include this third mode, we get 
the dashed line in Figure 4, which nearly sits on top of 
the solid line. In this case, the second excited mode 
gains less than 5% of the total population. 
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FIG. 6. This plot shows the fractional population in the |1) 
state for four different values of the detuning, obtained from 
a numerical solution of Eq ([!]). Starting from the top, the 
detuning is 6 = 0, 6 = 50 Hz, 5 = 100 Hz, and 5 = 200 Hz. 
The values of the other parameters are given in the text. 
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FIG. 7. This plot shows the period of modulation as a func- 
tion of detuning S. The dashed line corresponds to the numer- 
ical solution of the one-dimensional version of Eq. (|l]), while 
the solid line was obtained from a numerical solution of the 
two- mode model Eq. (|l8|). The Rabi frequency was fi = 700 
Hz. 



B. Dependence on detuning 

In Figure 6, we show how the behavior of the system 
depends on the detuning S. The Rabi frequency Q = 700 
Hz was held fixed for each plot while the detuning was 
varied from zero at the top 8 = to S — 200 Hz in the 
bottom plot. As predicted by the coupling parameter 



[3 — z 5Q/Ql s in the CGSD model, no coupling between 
motional states occurs if S = 0, and thus the Rabi oscil- 
lations experience no modulation. As S is increased the 
motional-state coupling becomes stronger and we expect 
the modulation frequency to increase. The amplitude of 
modulation also increases as the detuning is increased. 

We show the dependence of the period of modulation 
on detuning more explicitly in Figure 7. The dashed line 
is the numerical solution of the full problem given by 
Eq. (Q), while the solid line is the numerical solution of 
the two-mode model given by Eq. (pi 



C. Dependence on trap displacement 

In Figure 8, we show how the behavior of the sys- 
tem depends on the trap displacement zq. The Rabi 
frequency = 700 Hz and the detuning S = 100 Hz 
were held fixed, while the trap displacement was varied 
from zero zq = in the top plot to zq = 1 /im in the bot- 
tom plot. Again, the coupling parameter (3 predicts no 
modulation if = 0. As zq is increased, the frequency 
of modulation increases as the system is driven harder. 
However, for the large separation in the bottom plot, the 
modulation becomes highly irregular and the two-mode 
model most certainly breaks down. This behavior may 
be chaotic and warrants further investigation. 
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FIG. 8. This plot shows the fractional population in the 
|1) state for four different values of the trap displacement zo, 
obtained from a numerical solution of Eq (jl]) . Starting from 
the top, the displacement is zo — 0, za = 0.1 /im, zo = 0.4 /im, 
and zo = 1.0 /jm. The values of the other parameters are given 
in the text. 
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IV. CONCLUSIONS 

The gross features predicted by our model, such as 
double-peaked structure in the density distribution, and 
the presence of collapses and revivals in the relative pop- 
ulation dynamics, are supported by experimental obser- 
vation |25|] . Experiment-theory agreement on finer points 
is only fair. The theory tends to underestimate the con- 
trast ratio of the collapses, for instance. Moreover, to 
match the detuning trends shown in Fig. 6 and Fig. 7 
one needs to add by hand an unexplained overall detun- 
ing offset. This is most likely due to there actually being 
a spatial dependence of the bare Rabi frequency due to 
the influence of gravity on the untrapped intermediate 
state of the two-photon transition. To model the experi- 
mental situation in more detail one would have to include 
this effect as well as inelastic loss processes and finite- 
temperature effects neglected here. It may be also that 
treating the TOP trap potential as purely static may be 
an oversimplication. 

In this paper we have demonstrated the possibility 
for quantum state engineering of topological excitations 
through the interplay between the internal and spatial 
degrees of freedom in a Bose condensed gas. Due to the 
symmetry of the system we have analyzed, the excitation 
in our case was the odd-parity dipole mode. The intrigu- 
ing possibility of exciting modes with alternative sym- 
metries, such as a vortex mode [f3lH54|, would require a 
different trap geometry, but is a straight-forward exten- 
sion of the analysis presented here. Although we have 
focussed in this work on a particular parameter regime, 
the system is a rich one for study and exhibits complex 
and perhaps chaotic dynamics under strong excitation 
conditions. 
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TABLE I. This is a table showing the values used for the 
various physical parameters appearing in our calculations. 
The scattering lengths are taken from H. 
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